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Abstract 

The most recent crystal structure of the Fenna-Matthews-Olson (FMO) protein complex indi- 
cates that each subunit contains an additional eighth chromophore. It has been proposed that this 
extra site functions as a link between the chlorosome antenna complex and the remaining seven 
chromophores in FMO [Schmidt am Busch et al, J. Phys. Chem. Lett., 2, 93 (2011)]. Here, we 
investigate the implications of this scenario through numerical calculations with the generalized 
Bloch-Redfield (GBR) equation and the non-interacting blip approximation (NIB A). Three key 
insights into the population dynamics and energy transfer efficiency in FMO are provided. First, 
it is shown that the oscillations that are often observed in the population relaxation of the dimer 
composed of sites one and two may be completely suppressed in the eight site model. The pres- 
ence of the coherent oscillations is shown to depend upon the particular initial preparation of the 
dimer state. Secondly it is demonstrated that while the presence of the eighth chromophore does 
not cause a dramatic change in the energy transfer efficiency, it does however lead to a dominant 
energy transfer pathway which can be characterized by an effective three site system arranged in 
an equally-spaced downhill configuration. Such a configuration leads to an optimal value of the 
site energy of the eighth chromophore which is shown to be near to its suggested value. Finally 
we confirm that the energy transfer process in the eight site FMO complex remains efficient and 
robust. The optimal values of the bath parameters are computed and shown to be closer to the 
experimentally fitted values than those calculated previously for the seven site system. 
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I. INTRODUCTION 



The Fenna-Matthews-Olson (FMO) protein is one of the simplest and most well-studied 
light harvesting systems. The protein complex exists as a trimer of three identical subunits 
whose function is to link the chlorosome antenna complex, where light-harvesting takes place, 
with the reaction center, where charge separation occurs. FMO is also one of the earliest 
light harvesting systems for which a high resolution crystal structure has been obtained. 1 
This early crystal structure indicated that each of the three FMO subunits contains seven 
Bacteriochlorophyll (Bchl) chromophores which serve as the primary energy transfer path- 
way between the chlorosome and the reaction center.-"- Recently, however, a more careful 
crystallographic analysis of of FMO has been performed which demonstrates that the in- 
dividual subunits contain eight Bchls, not seven.-^ The eighth Bchl resides on the surface 
of the protein complex and it has been suggested that this additional chromophore is often 
lost during sample preparation. 

From an energy transfer perspective, the presence of an additional chromophore may 
challenge current understanding of how exciton transfer occurs in FMO. For example, in 
many previous studies on the seven Bchl complex, it is thought that two nearly independent 
energy transfer pathways exist. 8 " 12 Sites one and six are approximately equidistant from the 
antenna complex and both are assumed to be possible locations for accepting the excitation 
from the chlorosome. From there, the energy is subsequently funneled either from site one 
to two (pathway 1) or from site six to sites seven, five, and four (pathway 2). The terminal 
point through either route is site three where the exciton is then transferred to the reaction 
center (see Fig. |2fa)). The couplings within each of the pathways are much larger than the 
couplings between the two which implies that the two routes are nearly independent .—^^ In 
the second paper of this series, the two pathways and the respective probability of traversing 
each have been quantified using a flux analysis.— 

However, this two pathway picture is not entirely consistent with the recent experimental 
data. The new crystal structure indicates that the eighth chromophore resides roughly mid- 



way between the 



argued in Ref. 



Daseplate and the Bchl at site one.-^ 1 ^ Additionally, Renger and coworkers 



14 



that the eighth Bchl provides the most efficient path for exciton injection 
into FMO as a result of its position and orientation with respect to the chlorosome. If this 
is correct and site eight serves as the primary acceptor of excitation energy from the chloro- 
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some, then a preferential energy transfer route emerges through pathway 1. Due to the weak 
inter-pathway couplings, the secondary channel involving the remaining four chromophores 
in pathway 2 is largely bypassed in this scheme. This observation may have a significant 
impact on the efficiency and robustness of the energy transfer process. The main objective 
of this work is to address this issue by exploring how the dynamics and the energy transfer 
efficiency in FMO are affected by the presence of the eighth site and a realistic environment. 

The first major conclusion of the present study is related to the population dynamics 
in FMO. In many of the the previous studies of the seven site model of FMO, the popula- 
tion relaxation dynamics are modeled with site one or site six initially populated.- - — Under 
either of these initial conditions, pronounced oscillations in the short-time dynamics are 
observed. However when site eight is initially excited,— then the oscillations in the popula- 
tions are completely suppressed. This lack of oscillations has been independently observed 



in the dynamics recently reported in Refs. Il4j-|l6l Here, we provide a simple explanation 



for this behavior. The eighth Bchl maintains a large energy gap with the other seven sites 
in FMO in order to facilitate efficient directed energy transport. However, it is also rather 
weakly coupled to the remaining Bchls. This leads to a slow incoherent decay of the initial 
population at site eight, and hence a broad distribution of initial conditions at the dimer. 
The consequence of this result is that the population oscillations generally observed between 
sites one and two are completely suppressed which illustrates the importance of the initial 
conditions on the dynamics of the dimer. 

The second key result of this work demonstrates that if the eighth Bchl is the primary 
acceptor of excitation energy from the chlorosome as recently proposed then a primary energy 
transfer pathway in FMO does indeed emerge.— Note that this situation is substantially 
different from the previous interpretations of the energy flow in the seven site models where 
two independent pathways are generally assumed to exist. The extent of this effect and 
its impact on the energy transfer efficiency is quantified by introducing reduced models of 
FMO that consist of only a subset of the sites in the full system. It is demonstrated that 
sites eight, one, two and three which constitute pathway 1 provide the largest contribution 
to the dynamics of the full system. The remaining four sites of pathway 2 are seen to play 
a relatively small role. Despite the fact that only a single pathway dominates the energy 
transfer process, we also show that the presence of the extra Bchl does not significantly 
impact the efficiency or robustness of FMO. The eight site model leads to only a slight 
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increase in the transfer time as compared with the seven site system, and thus maintains 
the same high efficiency as observed in previous studies of FMO. 

Based upon energetic arguments, it has been suggested that the presence of the eighth 
Bchl leads to optimal energy transfer in FMO.— That is, its location near the chlorosome 
allows for a large coupling to the antenna complex as well as substantial overlap of the 
absorption spectrum of the eighth Bchl with the fluorescence spectrum of the chlorosome. 
These factors result in efficient transfer of the excitation energy into FMO while simultane- 
ously allowing the eighth chromophore to maintain a large energy gap with the remaining 
Bchls and hence a favorable energy transport landscape. These features implicitly suggest 
that there should be an optimal value of the site energy of the eighth Bchl. Here it is demon- 
strated that this observation is correct. However, in this case, the behavior is independent 
of the presence of the chlorosome, and can be understood by considering a further reduction 
pathway 1 to only three sites. The result of this procedure is a downhill configuration of 
three equally spaced sites (see Fig. [2(c) ) which is known to allow for highly efficient energy 
transfer.— 

Recently, several studies have shown that the environment does not have an entirely 
destructive role in the energy transport properties of excitonic systems.— d!d§r2H Instead, the 
environment can serve to enhance both the efficiency and robustness of the energy transfer 
process. Optimal values have been shown to exist for the temperature as well as other 
bath parameters which maximize the energy transfer efficiency in several light harvesting 
systems. Moreover, the experimentally fitted model parameters for FMO are near optimal 
in many cases. An extensive search for the optimal environmental parameters has been 
recently presented in Refs. Ill and 2l|. In addition to the above findings, we also explore 
the effect of the eighth Bchl on the environmentally assisted energy transport properties 
in FMO. It is found that the optimal values of the bath parameters are similar to those 
found in the seven site model, but are closer to the experimental values in general. As has 
been observed before, the energy transfer efficiency is relatively stable over a broad range 
parameters illustrating the robustness of the network.— 

In the next section, we present the average trapping time formalism which is used in 
the remainder of the discussion as a measure of the energy transfer efficiency.— 1 ^ This is 
followed by a brief outline of the generalized Bloch-Redfield (GBR) approach and the model 
Hamiltonian for the eight site FMO complex used in the numerical calculations. The results 
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for the population dynamics and the development of the reduced models for FMO are then 
presented in Sec. II III This is followed by calculations of the trapping time as a function of 
the site energies and bath parameters. There it is demonstrated that optimal values exist for 
many of these parameters, and additionally that the experimentally fitted values for FMO 
are near-optimal. 



II. METHODS 



A. Average Trapping time 



The formalism for calculating the averagin g tr apping time in light harvesting systems 



17 and 
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Here we provide only the salient 



has been presented in detail previously in Refs. 
results. The total system is characterized by a discrete iV-site system Hamiltonian, H s and its 
interaction with the environment, H s b. Each site of the system is coupled to an independent 
bath of harmonic oscillators with the respective Hamiltonians, H b — ~ ^\ [p 2 - + tOjX 2 ) . The 
total Hamiltonian is then given by, 

N N N 



n n^m n 



n\ 



H, 



(n) 



Jj 1 ) (n) 



where denotes the coupling coefficient of site n to the j-th mode of its associated 
bath. The values of the site energies, e n , and couplings constants, V^ m , are specified below 
in Sec. [HDl 

The time evolution of the reduced density matrix of this system can be conveniently 
described in the Liouville representation as, 

dp® 



dt 



-L tot p(t) 

-(L B + Ltr&p + -^decay + -^sb) p(t) j 



(2) 



where L s p = i/h [H s , p] describes the coherent evolution under the bare system Hamiltonian 
H s . In light harvesting systems, the energy flows irreversibly to the reaction center which is 
modeled here through the trapping operator [L trap ] nni)TlTO = (k tm + k tn ) /2, where k tn denotes 
the trapping rate at site n. The energy transfer in FMO exhibits almost unit quantum yield. 
As a result, the decay rate of the excitation at any site to the ground state, k d , is expected 
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to be much smaller than the trapping rate, k t ^> kj. This allows for the simplification 

L decay 0. 



B. Generalized Bloch-Redfield Equation 

It remains to account for the Liouville operator describing the system-bath coupling L s b 
in Eq. El For a harmonic bath linearly coupled to the system, the time correlation function 
of the bath coupling operators is given by the standard relation^^ 

'h/3u 



C{t) = ^J duj J (u) (coth 



cosiut) — i sin(wt) , (3) 



where j3 = l/k^T and J(u) = §X/j ^"^(^ ~~ u j) * s the spectral density of the bath. For 
simplicity, we assume the spectral density is the same for each of the independent baths and 
given by the Drude form 

J{u) = 2^-^—^ , (4) 

where A is the bath reorganization energy and u c is the Debye cutoff frequency. For this spe- 
cial choice, the correlation function may be expanded in terms of the Matsubara frequencies, 
v i - Tpi as— >— >— 



\ j=l J / j=l J 

oo 

= j> ie -*, (5) 

j=0 

which defines the complex expansion coefficients otj with the condition u = u c . 

The dynamics in FMO have been computed using a variety of methods ranging in both 
accuracy and cost.-*^^^ 1 ^ 1 ^ Here, we choose the approximate generalized Bloch-Redfield 
(GBR) method which follows from a second order cumulant expansion in the system-bath 
interaction. It provides an accurate, but computationally friendly approach for the propaga- 
tion of the density matrix over much of the physically relevant parameter space.— ^ Due to 
the decomposition of the bath autocorrelation function in Eq. the system-bath interaction 
may be accounted for through the introduction of auxiliary fields. The GBR equation of 
motion for the reduced density matrix is then given by 

= - (L sys + L trap ) p(t) - iJ2 E 9nj{t)] • (6) 

n j=0 



The coupling of each Bchl to an independent bath leads to the additional sum over the N 
sites where the system operator A n = \n) (n\ and g n j denotes the jth-auxiliary field coupled 
to site n. The auxiliary variables are subject to the initial conditions g n j{0) — and obey 
the equations of motion, 

d9n ^ = - (L sys + L trap + Vj ) g n>j (t) - iRe{ aj ) [A n , p{t)) + Im(a i ) [A n , p(t)} + , (7) 

where the plus subscript denotes anti-commutation. 



C. Trapping Time 

The mean residence time at site n is by definition 

/■oo 

( Tn ) = / dtp nn (t) , (8) 
Jo 

where p nn denotes the population at site n. The average trapping time is then given simply 
as the sum of the residence times at each of the N sites, (t) = J2n=i( T n)- Invoking the 
steady state solution of Eq. El L tot (t) = p(0), then the average trapping time is given by the 
compact expression, 

(t) = Tr (0(0)) , (9) 
where the trace is taken over the site populations of the reduced density matrix. 



D. Eight Site FMO model 

The Hamiltonian for FMO is constructed from the crystal structure recently deposited in 
the protein data bank (pdb code: 3eoj).- The site energies are taken from those computed in 



Ref. 



l4 and the coupling element between sites n and m is calculated from the dipole-dipole 



approximation, 

( d n • d m (d n • r nm ) (d m • r nm )\ 

Vnm = O I - ^ 6 ■ -5 I . {10) 

Additional details and the explicit system Hamiltonian are given in the Appendix. Aside 
from the eighth site, the most significant difference between the present model Hamiltonian 
and the model previously derived by Cho et al.— is in the energy difference between sites one 
and two. In the current case, the energy transfer through pathway 1 is entirely energetically 



favorable whereas a barrier is present between site one and two in the model of Ref. 
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Unless otherwise stated, the bath is characterized by the experimentally fitted values 
for the reorganization energy of 35 cm -1 and Debye frequency uj" 1 = 50 fs (105 cm" 1 ).^ 
Additionally the temperature is 300 K and the trap is located at site three with a trapping 
rate of k t = 1 ps _1 . 



III. NUMERICAL RESULTS 



A. Population Dynamics and suppression of the oscillations 

The time evolution of the populations in the eight site model of FMO calculated 
from Eq. [2] using the GBR is shown in Fig. Ufa) and (b) for the initial population located 
at site one and site eight, respectively. The bath parameters are taken at their fitted values 
specified above and the trap at site three is not included. The most striking difference seen 
between Fig. []Ja) and (b) is the absence of the oscillations in the populations of the dimer 
when site eight is initially excited. Other recent studies of the eight site model of FMO 
have also observed a similar lack of oscillations.— ~— The origin of this effect may be traced 
to the initial conditions at the dimer. The energy difference between site eight and the 
remaining sites is much larger than any of its respective couplings. This leads to the rather 
slow incoherent exponential relaxation of the population of site eight seen in Fig. QJb). The 
resulting initial conditions at sites one and two are then given by a corresponding incoherent 
distribution. It is this dephasing that suppresses the oscillations generally observed in the 
dynamics of the dimer. 

By applying the non-interacting blip approximation (NIBA) to the spin-boson model, 
Pachon and Brumer established that a necessary condition for the presence of the oscilla- 
tions in the dimer is an effective low temperature.™ The results of Fig. [T] demonstrate that 
the initial conditions impose an additional constraint on the observation of population os- 
cillations. Note that there are a variety of other initial preparations -such as starting from 
an eigenstate of the total system or exciting the system with incoherent light 2 -- which will 
also suppress the oscillations in the dimer. 

In order to analyze the influence of the initial conditions in more detail, the dynamics 
of the dimer calculated using the NIBA are presented in Fig. [TJc) and (d). The popula- 



tion dynamics described in Ref. 
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may be formulated as an equivalent generalized master 
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equatior 28 '" 



or,M) -' v 



f dt' J2 K nm(t-t')P m (t') , (11) 
J( > m=l 



dt 

where P n (t) denotes the population of site n at time t. The elements of the time- dependent 
transition matrix are constructed in the standard fashion 31 

K nm (t) = (1 - 5 nm ) W mn (t) - S nm W„fc(f) (12) 

fe 

where 5 nm is the Kronecker delta function and the individual rate kernels are given by the 
NIBA 

W nm (t) = 2 V* m e ^n-e m -2X)t-2C(t) _ (13) 

As defined previously, the coupling between site n and m is denoted by V nm , e n is the 
energy of site n, A is the reorganization energy, and C(t) is the bath correlation function 
given in Eq. |5j The results shown in Fig. [TJc) are calculated from Eq. [11] with the initial 
population located at site one and are seen to capture the key features of the full GBR 
dynamics shown in Fig. HJa). The decay is accounted for by setting the transfer elements 
K n i and K n2 to zero for all sites n > 2 which allows for population transfer from the dimer to 
the remaining Bchls, but prevents any back-transfer. Effectively this results in the addition 
of traps at sites one and two and thus leads to the population decay of the dimer. Without 
this decay, the two-site dynamics reproduce those of Ref. 
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As has been alluded to previously, the lack of oscillations in the dimer when the initial 
population is located at site eight may be explained by creating a distribution of initial 
conditions at site one. All of the population is initially located at site one in Fig. Ufa), 
whereas in Fig. [T](b), the corresponding initial conditions are given by the population that 
steadily flows from site eight. This slow incoherent decay of the population at site eight may 
be adequately fit to the single exponential Ps{t) = exp(— jt) where 7 = 3 ps _1 . Assuming 
that site eight decays only into site one, then the population of the latter is given by 1 — P§(t), 
and the corresponding transition rate is Wsi(t) = 7exp(— jt). The influence of site eight 
on the oscillatory behavior of the dimer can then be captured by convoluting the dynamics 
given in Fig. [T](c) (denoted by P n (t)) with this initial condition,— 

P n (t)= f dt' P n (t-t')W 81 (t') . (14) 
Jo 

The result of this procedure is shown in Fig. QJd). As is evident, the oscillatory behavior has 
completely disappeared. For large 7, the transition rate becomes a delta function and the 
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dynamics of Fig. [Tfc) are recovered. However, oscillations in the populations of both states 
can be observed only if the decay rate out of site eight is increased fivefold. The presence 
or absence of the trap simply effects the long time decay of the dimer populations and is 
irrelevant for the short time oscillatory behavior. These results demonstrate the importance 
of the initial preparation of the populations on the oscillations in the dynamics. 

It should be noted that while the NIBA calculations lead to qualitatively similar popula- 
tion dynamics as those given by the GBR, neither of the two approaches are exact. In FMO 
and other light harvesting systems, many of the model parameters are of the same order of 
magnitude. For instance, the couplings, V nm , and energy differences, e n — e m , as well as the 
reorganization energy, A, and thermal energy, are all of comparable magnitude. As a 
result, methods based upon second-order perturbation theory are, in general, insufficient to 
quantitatively describe the dynamics. A systematic procedure for computing higher-order 
contributions to the NIBA rates has been derived and recently implemented.— 1 ^ This leads 



to non-negligible corrections to the dynamics in t 



re spin-boson model, FMO and LH2. 
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25l capture the qualitative features that 



Thus while the results in Fig. [T] and those of Ref. 
are necessary to analyze the energy transfer behavior, there will be quantitative corrections 
from higher order terms. 



B. Pathway analysis and the ladder configuration 

Regardless of the presence or absence of oscillations, it is readily seen 
from Fig. (T](a) and (b) that the population primarily flows through pathway 1. Among 
the sites in pathway 2, only site four ever accumulates more than ten percent of the popula- 
tion. Particularly for times less than 500 fs, the sites from pathway 2 are scarcely populated. 
Similar behavior of the population dynamics has been seen in many other simulations of the 
seven site model for FMO.-^— These observations lead to the first reduced model for 
FMO which consists of only the four sites in pathway one shown in Fig. [2](b). As demon- 
strated below, this model is able to accurately capture the key features of the energy transfer 
process. One may proceed further by noting that sites one and two are coupled more strongly 
than the energy difference between the two. Additionally, sites eight and three are widely 
separated from either site in the dimer. The couplings between the distant sites and dimer 
are also rather weak (see values in the model Hamiltonian in Eq. [15] and Fig. (21(b)). As 



11 



a result, there can be rapid coherent energy transfer between sites one and two, whereas 
the transfer to sites eight or three will be comparatively slow. Therefore, when the initial 
population is located at site eight we may simply assume that these two sites of the dimer 
behave as one effective site with an energy that is the average of the two (270 cm -1 ). A 



simi 
Ref. 



ar "mean state" idea for developing this type of reduced model has been suggested in 
34 which explores the behavior of a dimer embedded in the PC645 photosynthetic net- 



work. Furthermore, the coupling between site eight and the terminal site, three, is negligibly 
small. This leads to a three site model for FMO where the couplings are determined by the 
nearest-neighbor values as shown in Fig. |2(c). 

C. Site Energy of Bchl 8 

Fig. [3] displays the average trapping time calculated as a function of the site energy of Bchl 
eight. It contains two additional key findings of this work. The first is that the trapping time 
behavior, and hence the efficiency, seen in the eight site model of FMO is largely governed 
by the sites in pathway 1. The second feature is that an optimal value exists for the site 
energy of Bchl 8, and moreover, the optimum is near the fitted value determined in Ref. Il4 . 
The source of the latter is the highly efficient ladder configuration shown in Fig. [2(c). 

The main portion of Fig. [3] displays the average trapping time calculated with the full 
Hamiltonian given in Eq. [TB] along with the corresponding results for the four site model 
of FMO (see Fig. [2(a) and (b)). The bath parameters are taken at their experimentally 
fitted values with the temperature of 300 K. The inset of Fig. [3] contains the results from 
the three sites model shown in Fig. [2(c). For this case, the exact trapping time calculated 
from the hierarchical equation of motion method is also presented, as well as the results of 
a Forster theory calculation. As can be seen, both the GBR and the Forster calculations 
predict optimal values of the site energy that semi- quantitatively capture the behavior of 
the exact hierarchical results. 

These results demonstrate that the energy transport is dominated by a subset of the sites 
in FMO and furthermore that the mechanism is correctly described by Forster theory. The 
four site model correctly describes the qualitative features seen in the full system and it 
accounts for the majority of the trapping time. Of the remaining sites in Fig. [TJ site four 
was seen to have the largest impact on the population relaxation. Calculations that consist 
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of pathway 1 plus site four are seen to capture almost all of the behavior seen in the full 
eight site system. 

In Ref. uM, it was noted that the site energy of Bchl eight maximizes the overlap with 
the chlorosome emission spectrum while simultaneously maintaining a large energy gap with 
the remaining seven core Bchls. This indicates that there should be an optimal value of the 
site energy of Bchl eight. In addition to the observation that the trapping time behavior 
may be captured by a simplified model of FMO, there is another interesting feature seen 
in Fig. |3j The trapping time displays a minimum as a function of the energy of the eighth 
site for all of the constructed models. Moreover, for the eight site model, the optimal value 
of the site energy is rather close to the fitted value of 505 cm" 1 .— Increasing the energy 
difference between site eight and site one decreases the back-transfer rate, but also decreases 
the spectral overlap between the two. The position of the optimal value is no coincidence. 
The three-site model in Fig. [2](c) readily demonstrates that this particular choice for the site 
energy leads to a downhill configuration of three sites that are approximately equally-spaced 
and equally-coupled which is known to be very efficient.— The qualitative behavior of the 
trapping time in the full, complicated eight site system becomes obvious with the aid of the 
reduced models. Note also that the average trapping time varies little over a large range of 
values of the energy of site eight indicating the robustness of the energy transfer process. 
Below it is demonstrated that many of the other fitted parameters are near optimal as well. 



D. Optimal Bath Parameters 

The average trapping time calculated as a function of the reorganization energy is shown 
in Fig. HJ^a) and (b). Fig. H^a) varies the reorganization of all eight sites simultaneously 
(in this model all chromophores are assumed to have identical environments) whereas (b) 
varies only that of site eight while keeping all of the others at the fitted value of 35 cm" 1 . 
In order to demonstrate that the presence of the additional chromophore does not lead to a 
dramatic increase in the trapping time, two different initial conditions are taken with either 
site one or site eight initially excited. As expected, a slightly larger trapping time is observed 
with the initial population at site eight due to the larger distance to the trap. However, 
the difference between the two scenarios is not substantial. Additionally, the qualitative 
behavior of the two initial conditions is quite similar and both lead to an optimal value of 
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the reorganization energy that is close to the experimentally fitted value of 35 cm -1 . It has 
been proposed from recent numerical simulations that the reorganization energy in FMO 
should be approximately twice as large as the experimentally fitted value used here.— The 
optimal values of A in Fig. 11(a) (55 cm -1 ) and (b) (40 cm -1 ) are consistent with a somewhat 
larger value of the reorganization energy. The mean trapping time is more sensitive to 
variations of the reorganization energy than was observed for the site energy in Fig. [3j but 
there is still a large range of A where the trapping time is near optimal. 

Finally, Fig. H^c) and (d) display the results for the average trapping time calculated as 
a function of the Debye frequency and as a function of the temperature. As in Fig. IH(a) two 
initial conditions are shown with the initial population at either site one or site eight. Again, 
the average trapping time from site one is always faster than for site eight. Nevertheless, 
the two initial conditions display qualitatively similar behavior for both the temperature 
and Debye frequency. Additionally an optimal value of the Debye frequency is observed 
that is close to the experimentally fitted value of 105 cm -1 . These results for the average 
trapping time as a function of the bath parameters are similar to those observed previously 
for the seven site model of FMO.— There is one notable difference however. In all cases, 
the Hamiltonian for the eight site FMO model recently proposed by Renger et at^ leads to 
optimal conditions that are closer to the experimentally fitted values than those calculated 



previously^ using the Hamiltonian of Ref. 
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IV. CONCLUSIONS 



The FMO protein serves as one of the model light harvesting systems, and the qualitative 
features of the energy transfer process have been understood for some time. However, recent 
experimental evidence has shown that many of the previously developed theoretical models 
are not entirely complete. An additional chromophore is present in each subunit of FMO 
that resides between the chlorosome and site one. In this work we have shown that the 
presence of the eighth site does not significantly alter the previous conclusions that have 
been reached with regards to environmentally-assisted exciton transport. 11 Optimal values 
exist for many of the bath parameters and, moreover, the optimal conditions are generally 
closer to their respective experimentally fitted values than in the seven site FMO models. 
Additionally, the dependence of the average trapping time with respect to variations of the 
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bath parameters is rather weak illustrating the overall robustness of the energy transfer 
process. 

However, the presence of the eight Bchl may necessitate a reassessment of our under- 
standing of the energy transport process in FMO. Given that site eight is the primary entry 
point for the exciton into FMO, then Fig. [T] clearly exhibits a complete suppression of the 
population oscillations that are generally observed in the seven site models of FMO. That 
is, the coherent population oscillations observed in previous studies depend upon the spe- 
cial choice of initial conditions. Here we have shown that the origin of the suppression is 
the slow decay of the initial population at site eight which leads to an incoherent distribu- 
tion of initial conditions at the dimer. In the physical setting there will be an additional 
source of dephasing due to the extra step from the chlorosome to FMO. This will broaden 
the distribution of initial conditions even more and further suppress the oscillations in the 
dimer. 

An additional feature of the eight site model that is markedly different from the previous 
seven site configuration is observed in the population dynamics shown in Fig. [1] and the 
average trapping times displayed in Fig. |3j These results demonstrate that the energy flow 
in the eight site model is dominated by a subset of the chromophores, whereas it has been 
previously assumed that two independent pathways involving all of the Bchls are available 
for the energy transfer process. The qualitative features of the transport in the eight site 
model are largely determined by the dynamics of pathway 1. Sites four, five, six and seven 
provide a rather small contribution to the overall efficiency in this case. The agreement 
between the results for the full eight site system and the reduced four- and three-site models 
shown in Fig. [2] provide further support to this claim. Nevertheless, the eight site model and 
the seven site model display similar energy transport efficiencies. The origin of this behavior 
in the former is evident from Fig. [2](c) which shows that the eighth Bchl forms an optimal 
downhill ladder configuration with the dimer and site three. This result demonstrate the 
usefulness of the reduced models in providing an intuitive explanation of many of the key 
features present in the numerical results. 
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VI. APPENDIX 



Following the prescriptions used previously for constructing the dipole-dipole 
interactions,— 1 ^ 1 ^ the unit vectors, d n , in Eq. [10] point along the axis connecting the iV& 
and Nd atoms of the n-th Bchl and r nm is the vector connecting the Mg atoms of Bchl n 
and m. Setting the constant C = 155000 cm _1 A 3 leads to an effective dipole strength of 30 
D 2 . With these specifications, the system Hamiltonian (in cm -1 ) for the eight site model is 



H, 



FMO 



( 310.0 


-97.9 


5.5 


-5.8 


6.7 


-12.1 


-10.3 


37.5 \ 


-97.9 


230.0 


30.1 


7.3 


2.0 


11.5 


4.8 


7.9 


5.5 


30.1 


0.0 


-58.8 


-1.5 


-9.6 


4.7 


1.5 


-5.8 


7.3 


-58.8 


180.0 


-64.9 


-17.4 


-64.4 


-1.7 


6.7 


2.0 


-1.5 


-64.9 


405.0 


89.0 


-6.4 


4.5 


-12.1 


11.5 


-9.6 


-17.4 


89.0 


320.0 


31.7 


-9.7 


-10.3 


4.8 


4.7 


-64.4 


-6.4 


31.7 


270.0 


-11.4 


v 37.5 


7.9 


1.5 


-1.7 


4.5 


-9.7 


-11.4 


505.0 j 



(15) 



where the zero of energy is 12195 cm 1 . Note that there is an error in the sign of the coupling 



between sites one and two in the table provided in Ref. [l4j. Aside from this, these values 



reproduce all of the couplings listed therein to within 3 cm 1 . 
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Figures 




Time (ps) 

FIG. 1: Site populations of in the eight site model of FMO with Bchl one (a) or eight (b) initially 
excited calculated with the GBR. The populations of the remaining sites five, six and seven are 
never larger than 10% and not shown. The site populations of the dimer calculated using the 
NIB A of Eq. 1111 calculated with site one initially excited (c) and from Eq. [TH (d). In all cases, the 
temperature is 300 K with a reorganization energy of 35 cm -1 and cutoff frequency of uj~ 1 = 50 
fs. 
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FIG. 3: The trapping time as a function of the site energy of site eight. The solid (black) line 
and dashed (blue) line in the main figure correspond to the results calculated from the full eight 
site Hamiltonian of Eq. [JS] and with only the four sites of pathway 1, respectively. The red dots 
correspond to the optimal site energies and the vertical dashed line indicates the fitted value of 



14. The inset contains results for the 



the site energy of Bchl 8 of 505 cm -1 determined in Ref. 
three site model calculated with the GBR (solid black line), Forster theory (dotted red line), and 
hierarchical equation of motion (dashed blue line). The remaining parameters are the same as 
in Fig. [TJ 
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FIG. 4: The trapping time as a function of the reorganization energy of all Bchls (a), and as a 
function of the reorganization of site eight only (b) while the remaining seven sites are fixed at the 
experimentally fitted value of A = 35 cm -1 . The trapping time as a function of the bath cutoff 
frequency and temperature are shown in figures (c) and (d) respectively. In all cases the trap rate 
at site three is 1 ps. The solid (black) and dashed (blue) lines correspond to initial excitation at 
site eight or site one, respectively. The red dots indicate the optimal trapping times and the dotted 
vertical lines correspond to the respective experimentally fitted values. The remaining parameters 
are the same as in Fig. [TJ 
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